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SUMMARY 


This final report consists of two reports concerning analytical tools. 

The first report Is entitled, ’'Probability Analysis of the Orbital Distribution of 
Events". It was included as an Appendix in AMA Report No. 73-44 (October 1973). 
The second report is entitled, "Particle Trajectories for a Repulsive Force Field" 
These two reports summarize the work done on Task III of the subject contract. 

The work on Task IV was included as a modification to the program re- 
ported in AMA Report No. 73-44 (Statistical Information Program, Progress Re- 
port, (Task II)), October 1973. 

The modification works fully automatically and points out TOF events. The 
modified program has been checked thoroughly and delivered to the Technical 
Officer. So far, two TOF events were discovered. 

The work on Task III and IV of the subject contract has thus been com- 


pleted and will be Invoiced. 
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PROBABILITY ANALYSIS OF THE ORBITAL 
DISTRIBUTION OF EVENTS 


1. Introduction. 

The fine structure of the dust flux density may be examined by noting 
asymmetries in the event rate. The analysis contained herein permits a dis- 
tinction between asymmetries caused by random fluctuations and systematic 
variations by calculating the probability of any particular asymmetry. The 
analysis is applied to an asymmetry observed in Pioneer 8 and 9 distributions. 

2. Analysis. 


Let the relative probability of an event occurring in the half orbit 
as the spacecraft is moving toward the solar apex be p, and in the other half 
orbit H be q, where p + q=l. 

A 

If n events are observed altogether, the probabilities of all possible 

combinations are given by the terms of the binomial expansion 

. , v n n /n\ n-1 /n\ n-2 2 

(p+q) =p +( il/ )P Q + ^Jp Q + ’** 

Thus the first term in the expansion gives the probability of all events occurring 
in H^, the second term the probability of all except one of the events occurring 
in H^, etc. Other probabilities may be obtained by summing appropriate terms 
in the expansion (1). For Instance, the probability of at least y of the events 
occurring in H (i.e., at most n -y events in H ) is given by the sum 

T “fV s »- k k A 

) U p “ • 

k=0 

It is apparent that for large values of n, the computation of the numerous terms 

IrjTTAlTrorl In /9 \ ma tt Kp o v't t*q m aNt loKnn’rtnc UAnrmmr fnr n nn » *1 on oHoamq+a 
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approximation to the binomial distribution is given by:: 


O 


k n-k 

P q 


(See Reference 1); 


«/2ffnpq 


-(k -d p) /2npq 


( 1 ) 


( 2 ) 
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i.e. , the binomial distribution may be approximated by a Gaussian distribution 
of mean n p and standard deviation vnpq . 

Instead of summing the terms involved, in an expression such as (2), we 
may use the corresponding extensively tabulated integrals of the normal curve. 


jlication to Actual Orbital Distribution. 


The actual distribution for Pioneer 8 events was found to be: 


Total number of events analyzed: 230 

Events in H^, (moving toward s.a. ) 139 


Events in H . 91 

A 

Corresponding numbers for Pioneer 9 

Total number of events analyzed: 111 

Events in 73 

Events in H. 38 

A 

Combined totals: 


Total number 341 

in H t 212 

in H 129 

A 

As a working hypothesis it is assumed p = q = .5 (i.e. , the probability of 
an event occurring is the same for each half orbit). 

For Pioneer 8. The expected number in each half orbit is then np = 115 
and the deviation from this expected number is 139 - 115 ~ 24, thus for this case 

a = ^230 *2 * 2 = 7,5S * Consequently, we have a deviation of 3.17 a, corres- 

ponding to a probability of 0.16%. 
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For Pioneer 9 the deviation from the mean is 17.5 and the standard deviation 
is 5.27 resulting in a deviation of 3.32 a. This corresponds to a probability of 
0.08%. 

It is realized, of course, that the time spent in each of the half orbits is 

not the same. Pioneer 8 spends 578 days in H and 507 days in H corres- 

1 A 

ponding to a daily event rate of . 240/dayand . 179/day, respectively. If we in- 
crease the time in H by 71 days and hypothesize that the daily event rate holds 

this would indicate 13 additional events occurring in the half orbit H . 

A 

So, the new hypothetical distribution would be: 

(Pioneer 8) 

Total number of events 243 

Events in H 139 

Events in H . 104 

A 

np = 121.5 
d = 17.5 
O = 7.79 

Thus d = 2.24 0, corresponding to a probability of 2.5%. 

For Pioneer 9, the event rate in H was .113 and the number of days 

A 

spent in FL, and H was 355 and 335, respectively. 

T A 

Making a same adjustment as above the additional 20 days would lead 
to an additional 2 events in H . Thus, the new hypothetical distribution would be: 

A 

(Pioneer 9) 


Total number of events 

113 

Events in 

73 

Events in H. 

A 

40 


Expected number, np = 56.5. 
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With the deviation d = 16.5, and the standard deviation a~ 5.32, then 


d = 3.10a. This corresponds to a probability of 0.2%. 

The combined results are tabulated as: 

Hypothetical: 

Pioneer 8 and 9 

Total number of events : 356 

Events in H„, 212 

T 

Events in H . 144 

A 

Expected number of events in each half orbit: 178 

d = 212 - 178 = 34 
a = 0. 43 
d = 3.61 a, 

corresponding to a probability of 0.04%. 

4. Conclusions. 

It is seen that the probability for random variations to produce the 
observed asymmetry in the distribution of events is very small. Then, the 
observed asymmetry should be attributed to some systematic effect, such as 
the presence of interstellar grains. 

5. References . 

[l] "Probability, Random Variables, and Stochastic Processes", 
Athanaslos Papoulls (Polytechnic Institute of Brooklyn). 

McGraw-Hill Series in Systems Science. Copyright 1965. 
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Introduction 


Under simple assumptions the force on a particle due to solar radiation 
pressure Is directed along the particle's radius vector, from the sun, and Is In- 
versely proportional to Its distance from the sun. In addition the force Is assumed 
proportional to the cross sectional area of the particle. Thus the particle's 
acceleration due to radiation pressure Is given by: 


C 0 A 
m 


R_ 

3 


and the equation of motion of the particle subject to both solar radiation pressure 
and gravitational attraction Is given by: 

" R °0 A R „ R 

HT— T^ n> 1 

r r r 


where: 


d) 


R Is the radius vector of the particle 
r Its magnitude (*|r [ ) 

2 

fig the sun's gravitational constant (= K M Q ) 

Cq the force due to solar radiation pressure on a particle 
of unit cross sectional area at unit distance 


A 

m 

fi 


the cross sectional area of the particle 
particle mass 
j©_ A_ 


If the particle is assumed approximately spherical then 

2 


A _ir a‘ 


where 


thttS 


m 4 3. 

3 9 


3 

4ap 


a 

P 


Is the particle radius and 
Its density 


3C 


O 


Pa 

defining the constant k as : 


(mp 2 ) l/3 


( 2 ) 


( 3 ) 


( 4 ) 


1 



/9ff\l/3 _©_ 


*( 


16/ 


For large, dense particles j3« 1; consequently the normal Keplerlan 
orbit for such a particle would only be slightly modified. As the particle radius 
or density decrease 0 will increase and may, for sufficiently small, light 
particles, exceed 1. In this case the net force will be repulsive. 


Kepler's theory may be modified for these cases. The motion then takes 
place along the branch of the hyperbola which is convex toward the focus occupied 
by the sun. 

For situations such as these there will be a ’forbidden zone" about the 
sun into which these particles may not penetrate. The zone is of course de- 
pendent on the "velocity at infinity", the radiation pressure and the angular 
momentum of the particle. It is the purpose of this note to develop a sketch 
of this "forbidden zone". 


The value of the radiation pressure constant used here, has been chosen 
to be consistent with that given by Sllverberg. {Reference [l] ). 

The attractive and repulsive cases are separated by 0-1; since here the 
particle moves in a force free field. For a particle of given mass this gives rise 
to a ’limiting" density; or, for a particle of given density, to a ’limiting" mass or 
diameter. Particles with densities (or masses) smaller than their limiting values 
are repelled, others are attracted. 

-13 

As an example: for a particle of mass 8. 9 x 10 grams ( similar to that 
for event 20 from Pioneer 8) the limiting density of 0.953 gr/cc. 

For a particle with a (very likely) density (p) of 3 gr/cc, its limiting dla- 
-5 

meter would be 3.85 x 10 cm, and it would have a corresponding limiting mass 
-14 

of 8. 99 x 10 gr. 

For the curves plotted and presented herein, a density of 3 gr/cc and an ex 
cess speed (V w ) of 20 km/sec (corresponding to the velocity of the solar system) 
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has been assumed. (Of course, corresponding curves for other values of these 
parameters can readily be obtained). 

A word description of the trajectory '‘flown " by a test particle In the presence 
of the assumed solar mass --radiation Influence is given below. In addition, there 
Is a sketch included here to show several geometric properties of the trajectory. 

The flight paths of interest are those which describe hyperbolae in the 
vicinity of the sun. These particles are assumed to "arrive" in the solar sys- 
tem with an excess speed (V^) of 20 km/sec (the assumed value). If these 
particles approach the sun not along a direct radius, but at some offset distance 
(i.e. , with non-zero angular momentum), then they will describe hyperbolic arcs 
which can be predicted from the two-body model described above. On the sketch, 
below, r V Is the closest solar approach distance; c Is the off-set distance 
(measured from the solar apex radial), and V # is the approach excess velocity. 
Also, as shown on the figure, 0 is an angle locating the hyperbola's axis 
measured from the solar apex radial. Due to symmetry, these hyperbolae may 
be oriented symmetrically about the radial. Thus there could be mirrored 
Images of each path on this principal radius. 

Computational Sequence 

In determining geometric characteristics of these trajectories, and ob- 
taining a description of the FORBIDDEN ZONE, the following calculations pro- 
cedure has been utilized. 


(1) . A density of p * 3 gr/cc is assumed; knowing (the limit 

mass), choose values for the particle mass (m < m^^). 

-5 

(2) . With k known, (9. 32 * 10 ) and having chosen p then calculate 

a j8 corresponding to each m; l.e., from (4) 


(P m) 


3 




Sketch depicting a particle (P) moving from », with a velocity V B , offset at a 
distance c, which has a closest approach (q) to the sun (fi), as shown. 
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Note : Here, 0 * 1. 0, by definition. 

(3) . Having obtained 0, calculate the equivalent p (sun value) for the 

system : 

with u = 2.959 * 10 -4 (AU 3 /DA 2 ). 

© 

( 4 ) . Next, assuming V oo = 20 km/sec, calculate the parameter, B, where: 

B 2 * 

V 

oo 

And, (5). Knowing B define characteristics for the hyperbolic path; e.g. : 

e * tan" 1 ( - tan' 1 (|) 

and 

a s ^ + [l+ (^) 2] i/2 

oo 

= B + [l+B 2 ] l/2 . 

These data may be used to plot a boundary defining closest approaches to the sun, 
for various system constraints and/or conditions. However, the FORBIDDEN 
ZONE, as defined, Is not yet obtained. This region of "space" lies on the sun's 
side of the geometric envelope for the family of hyperbolae traced out by the 
assumed particles. In this regard, then, it Is necessary to examine the pro- 
blem (mathematically or graphlcaUy) and to define the envelope accordingly. 

Since the model used for the present investigation leads to analytical re- 
sults, the envelope curve can be described mathematically. For this description 
the parameter (A) is Introduced, wherein A (see above), and the equations 
for the trajectory are rewritten. Making use of the formal operations to deter- 
mine an envelope, it can be shown that the pair of parametric cartesian equations 
for this curve are: 
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(A ) 2 [l -/ 1+A* ] - y 2 + 2 (AC) 2 [l + yi + A 2 ] 

X * — ^ — — — 

<Ay)[ /l + A S - 2 ] 

^ , (2 + 3A 2 )[l + J\ 

y = ±e { 

[l*A* 

For these expressions the x-axis parallels the solar-apex radial; and Its positive 
direction is in the apex direction. The y-axls is normal to x; it may be assigned 
signs of ± arbitrarily. (Recall that c is the off-set distance — for — ex- 
pressed as a distance, measured in AU. Consequently, the coordinates (x, y) are 
stmllarlly dimensioned). 

Typical Results 

Table I shows: (1), a summary for the calculations outlined above. Here 
one will find tabulated, as a function of m, the corresponding values for /3 and 
JJ (see Eqs. (2) and (3)). 

Also, in the table, is a listing of closest approach distance * or 

each of the m values selected. The quantity q^.^ (in AU) denotes the closest 
distance that a particle of mass m can approach the sun under the conditions 
assumed. Necessarily, these distances correspond to an offset distance (c) of 
zero. 

Next, for a unit offset distance (i.e. , c - 1.0 AU), typical values of 
q(AU) and Q (deg) are tabulated. These data are illustrative of the results which 
may be expected from the calculations procedure outlined earlier. 

Lastly, there is a listing of c(AU) values which are required to attain a 
q 25 2.0 AU (for each value of m defined). Since q * 2 AU has been arbitrarily 
chosen as an upper bound of interest, here, then these data define a practical 
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upper limit for the present investigation. (Note : when the notation c < 0 is used 
herein it infers that these particles cannot have a q . as small as 2.0 (AU), for 
the conditions implied). 

Figure I is a graph of q . and c, as functions of m, for the information 

min 

described above. 

Figure II is a polar plot of (q (AU), 0°) for values of /x (or m) based on 
a fixed off-set distance (c) of 1 AU. Since the curves are symmetric about 0° = O, 
then c may be the off-set on either side of the radial (0°); this Is a typical result 
for any value of c (within the range of possible q distances allowed). 

Figure III shows q, 0 boundaries for several of the values of m used 

here. That is, each boundary corresponds to a chosen value of m( s m ). 

lim 

This describes a limit in q for the allowable range in c (AU), at a fixed p. It 

should be remembered that q = 2.0 AU represents the arbitrary upper limit for 

-15 

this quantity. Note that when m < 30 * 10 gr then q > 2. 0 AU; i.e. , the 

particle is subjected to a large enough radiation (back-pressure) so that it does 

not have sufficient kinetic energy to move in as close as 2.0 AU from the sun, 

even for a direct radial approach. (The direct solar approach describes the 

q noted earlier). It should be apparent that if V >20 km/sec then the 
min 00 

corresponding q -values should be smaller ; the particles would be able to move 
closer to the sun against radiation pressure. 

Figure TV graphically depicts the FORBIDDEN ZONE for this problem's 
conditions and constraints. Each curve, sketched on the figure, corresponds 
to a selected value of m (< m^.^). These arcs describe the geometric envelopes 
mentioned earlier; they separate regions, in the vicinity of the sun, into which 
the particles may move from those in which they cannot (the FORBIDDEN ZONE). 
Due to the obvious symmetry of this situation these arcs are mirrored in the 
solar (apex) radial line. 
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TABLE I. Calculated Data 


Assume 

m(gr) 

Compute 

0 

Compute 

M(AU.DA) 

Closest Approach* 

<AU» 

Representative Data 
(Based on e * 1AU) 

q <AU) 61° 

c(AU> 

(For q = 2AU) 

89. 9508 * 10" 15 

1.0 

0 


1.0 

o 

o 

a* 

2.0 

85 * 10" 15 

1.01905 

-5,63685 * 10" 6 

0. 0845 

1.043 

87.58° 

1.956 

80 * 10~ 15 

1.03985 

-1.17924 * 10 -5 

0.1768 

1.0923 

84. 95° 

1.91 

75 * XO -15 

1.06247 

-1.84834 * 10~ 5 

0.2771 

1.1481 

82.11° 

1.86 

-15 

70 * 10 

1.087182 

-2.57973 ♦ 10~ 5 

0. 3867 

1.2119 

79.05° 

1.80 

65 * 10" 15 

1.114373 

-3.38430 * 10“ 5 

0. 5074 

1,2854 

75.76° 

1.73 

60 * 10“ 15 

1.144506 

-4.27593 * 10" 5 

0.641 

1.3706 

72.23° 

1.648 

ee 

55 * 10 

1.178187 

-5.2755 * 10" 5 

0.7909 

1.4708 

68.42° 

1.55 

50 * 10" 15 

1.216219 

-6.39792 * 10~ 5 

0. 9592 

1.5887 

64.38° 

1.448 

45 * IQ' 15 

1.259692 

-7.68428 * 10~ 5 

1.152 

1.730 

60. 06° 

1.304 

40 * 10" 15 

1.310132 

-9.17681 * 10 -5 

1.376 

1.9017 

55.48° 

1.12 

35 * 10‘ 15 

1.36976 

-1.094132 * 10” 4 

1.640 

2.1135 

50. 64° 

0.85 

30 * 10~ 15 

1.44199 

-1.30784 * 10" 4 

1,9607 

2.3808 

45.57° 

0.294 

25 * 10" 15 

1.53234 

-1.57519 * 10‘ 4 

2.3615 

2.7281 

40. 26° 

<0 

20 * 10’ 15 

1.65066 

-1. 92531 * 10” 4 

2.8865 

3.2 

34. 72° 

<0 


♦Closest Approach would correspond to c = 0(AU). 





Figure I. Closest approach (q m ^ n ) (on a radial) and offset distance (c), for 
q - 2 AU, as a function of mass (gr). 
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30 20 iO O & 20 30 

Give*) 

Figure II. Description of closest solar approach fq, 6) as a function of particle 
mass (m) - or equivalently (/? or (i) - for a fixed offset distance (c) = 
1.0 AU. 
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Figure III. A plot of q and 0for various values of c (c not shown). Note: each 
curve corresponds to a value of m, as specified. 
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